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ABSTRACT 

The first galaxies in the Universe are built up where cold dark matter (CDM) 
forms large scale filamentary structure. Although the galaxies are expected to 
emit numerous Lya photons, they are surrounded by plentiful neutral hydro- 
gen with a typical optical depth for Lya of ~ 10 5 (H I halos) before the era of 
cosmological reionization. The H I halo almost follows the cosmological Hub- 
ble expansion with some anisotropic corrections around the galaxy because of 
the gravitational attraction by the underlying CDM filament. In this paper, we 
investigate the detectability of the Lya emissions from the first galaxies, exam- 
ining their dependence on viewing angles. Solving the Lya line transfer problem 
in an anisotropically expanding H I halo, we show that the escape probability 
from the H I halo is the largest in direction along the filament axis. If the Lyo; 
source is observed with a narrow-band filter, the difference of apparent Lya line 
luminosities among viewing angles can be a factor of > 40 at an extreme case. 
Furthermore, we evaluate the predicted physical features of the Lya sources and 
flux magnification by gravitational lensing effect due to clusters of galaxies along 
the filament. We conclude that, by using the next generation space telescopes 
like the JWST, the Lya emissions from the first galaxies whose CDM filament 
axes almost face to us can be detected with the S/N of > 10. 
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1. INTRODUCTION 

One of the most important questions in cosmology is how galaxies are formed and 
evolved in the context of the cold dark matter (CDM) Universe. To answer the question, 
it is essential to search for young galaxies at high-redshift systematically and to study their 
observational properties in detail. In the hierarchical structure formation scenario of the 
CDM Universe, the first bright objects and/or the first galaxies are predicted to be formed 
in about some hundreds Myrs after the Big Bang (i.e., redshift ~ 10) and to be embedded 
in the deep gravitational potential wells of the CDM (see e.g., Abel et al. 1998; Bertschinger 
1998; Bromm et al. 1999; Yoshida et al. 2003). They are expected to be luminous at the Lya 
line (Partridge & Peebles 1967) because the strong ionizing radiation from young massive 
stars in the galaxies should lead to prominent Lya emission through the recombination of 
hydrogen in their interstellar medium. 

In this decade, many galaxies which appear to be at early stages have been observed 
(e.g., Rhoads et al. 2000; Steidel et al. 2000; Ajiki et al. 2003; Shapley et al. 2003; Hayashino 
et al. 2004; Matsuda et al. 2004; Wang et al. 2005). Recently, the galaxies beyond redshift 6 
have been frequently detected by their prominent Lya emission with faint continuum (Lya 
emitters [LAEs]) by using ground-based large telescopes like Subaru (Kodaira et al. 2003; 
Nagao et al. 2004; Taniguchi et al. 2005; Shioya et al. 2005), Keck (Hu et al. 2002; Rhoads 
et al. 2004) and VLT (Cuby et al. 2003; Kurk et al. 2004; Tran et al. 2004). Moreover, 
surprisingly, the possible detection of a redshift ~ 10 galaxy is reported by Pello et al. 
(2004), although this is still controversial (Bremer et al. 2004; Weatherley et al. 2004). In 
the following decade, next generation space telescopes like the James Webb Space Telescope 
(JWST) are launched. The expected limiting 10cx flux of JWST for a point source reaches 
to nJy-level (~ 10~ 31 erg s _1 cm -2 Hz -1 ) in the wavelength range of ~ 1 — 5 jttm; this range 
corresponds to the redshifted Lya wavelength at the source redshift of 7 — 40. Thus, it seems 
to be in not-so-distant future that the Lya emission from the first galaxies and/or the LAEs 
beyond redshift of 10 are firmly detected. 

However, before the era of cosmological reionization, these first galaxies should be sur- 
rounded by plentiful neutral hydrogen in the intergalactic medium (IGM); we call this "H I 
halo" from now on. For completely neutral Universe, typical scattering optical depth for the 
Lya photons escaped from the galaxies at redshift about 10 is estimated to be ~ 10 5 by us- 
ing the standard cosmological parameters. One might be discouraged by the extremely high 
optical depth for the Lya photons, but it has been shown that most of them can escape from 
the H I halo around a Lya source at young Universe and travel freely toward the observer 
(Loeb & Rybicki 1999). The reasons are as follows: (1) Lya photons are not destroyed in 
the intergalactic space because free-electrons and dust as sources of the two-photon decay 
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and the ultraviolet continuum absorption, respectively, rarely exist in the IGM in such early 
Universe; (2) in the H I halo, which is assumed to expand following the pure cosmological 
Hubble flow, the Lya photons are scattered many times and diffused redward in wavelength; 
at last, cumulative scattering redshift grows large enough to escape from the H I halo. 

The hierarchical structure formation scenario of the CDM Universe indicates that the 
first galaxies are formed at high-density peak regions in the CDM distribution. In such 
regions, the CDM structures are predicted to be filamentary by numerical N-body simulations 
(e.g., Yoshida et al. 2003). The underlying CDM filament may affect the H I halos around 
them by gravitational attraction force. It is naturally expected that the gravity works 
anisotropically; it is the strongest in direction perpendicular to the axis of the CDM filament 
(r-direction), while it is the weakest along the axis direction (z-direction). This situation 
around one of the first galaxies is illustrated schematically in figure 1. Considering such 
situation, we assume an improved expansion law to the H I halo in our previous paper 
(Kobayashi & Kamaya 2004, KK04 in short) (spheroidal expansion law; see equation (1) 
below): the H I halo expands following a pure Hubble expansion in z-direction, while it is 
decelerated to some extent from the pure Hubble flow in r-direction. 

We have investigated the effects of the anisotropic expansion of the H I halo to the Lya 
line profiles and luminosities of the first galaxies in KK04. We have found the following two 
effects: (1) the peak of the Lya line is redshifted to longer wavelength; (2) the full width 
at half maximum (FWHM) of the line is broadened to wider than those of the isotropic 
expansion model to the H I halo. Since the FWHM corresponds to the Doppler width at the 
temperature of ~ 10 8 K, it has been predicted that no intrinsic Lya line features of the first 
galaxies appear. As deceleration against the Hubble expansion by the gravity of the filament 
becomes strong, these two effects result in less luminosity of the diffuse Lya emission line. 

Furthermore, as expected from the anisotropy in the expansion of the H I halo, the Lya 
line profile has dependence on the direction from which the Lya source galaxy is observed. 
Clarifying such dependence is essential to discuss the detectability of Lya emission lines from 
the first galaxies. The dependence can be studied by dividing the solid angle around the 
source into some regions and by investigating the Lya line profile in each region. However, we 
have not been able to examine the dependence in KK04 to keep the high statistical reliance of 
the numerical solution in a Monte Carlo simulation. This is because the following reason: in 
a Monte Carlo simulation, the convergence of the solution needs to be checked by increasing 
the number of test photons, but the total number of test photons has been restricted by 
the weak computational power to < 10 8 in KK04. We have improved the numerical code to 
solve the transfer problem with larger number of test photons than ~ 10 8 . Therefore, this is 
the time to revisit the Lya line transfer problem in the anisotropically expanding H I halo 
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with the improved numerical code. 

In this paper, we reexamine the Lya line transfer problem with 10 times larger number 
of test photons than that of KK04. In §2, followed by KK04, we summarize the formalism of 
the Lya line transfer problem in an expanding H I halo and discuss about the model of the 
H I halo. Our numerical results such as Lya line profiles, apparent Lya line luminosities, 
and typical angular sizes of the first galaxies in Lya emission observed from the z- and 
r-directions are shown in §3. Then, in §4, we give the Lya flux of the source with intrinsic 
Lya line luminosity of 10 43 ergs s" 1 , and discuss the expected signal-to-noise ratio (S/N) for 
the source in the case of the observation by the JWST. Finally, we conclude the paper in §5. 
We adopt the WMAP cosmological parameters in the whole calculations of this paper, that 
is, tt b = 0.044, tt M = 0.27, fi A = 0.73 and h = 0.71 (Spergel et al. 2003). 



2. FORMULATION AND CALCULATION 

As briefly described in introduction, we have formulated a Lya line transfer problem 
in KK04 according to Loeb & Rybicki (1999), setting a simple assumption to the expansion 
law of H I halo, that is, spheroidal expansion law: 

v = H s [e*(x 2 + y 2 )+z*] 1/2 , (1) 

where v is recession velocity, H s is the Hubble expansion rate at the Lya photon source 
redshift of z s , and e is the asphericity or the decreasing rate of the recession velocity in 
r-direction (x — y plane) to the pure Hubble expansion. In other words, this parameter e 
represents the relative intensity of the gravitational attraction force by the underlying CDM 
filament to the cosmological Hubble expansion. We choose the parameter of e as a free 
parameter from 1.0 (corresponds to the pure Hubble expansion in all directions) to 0.3 with 
a step of 0.1. As noted in KK04, e is assumed to be a single value globally for simplicity, 
although it should have different values depending on the distance from the CDM filament 
and approach to unity as leaving the galaxy far behind. This assumption can be hold because 
of the following reason. The properties of escaping Lya photons are determined essentially at 
a position where the optical depth for the Lya photons decreases to about unity and where 
the physical distance from the source galaxy is ~ r*. The characteristic proper distance of 
is a function of e and gradually increases from 1.0 to 2.9 Mpc as e decreases from 1.0 to 
0.3. The gravity of the underlying CDM filament can affect this region even at the early 
Universe (Bond et al. 1996; Colberg et al. 2004). 

Let / = I{v, fi, r) be the specific intensity of a Lya line (in photons cm -2 s -1 sr" 1 Hz" 1 ), 
where v is comoving frequency, r is radial vector from the source galaxy, and O is direc- 
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tion cosine vector. The comoving transfer equation for a Lja line with the assumption of 
isotropic coherent scattering is given by 

dl 

n-VI + a^(e,w)— = xAJ-I) + S- (2) 

Here, v is the redshift of the frequency v a — f p hoton where v a is the resonant frequency of Lya 
and ^photon is the photon frequency; \v is the Lya scattering opacity at the frequency redshift 
v (the analytical form is represented in KK04); J = (1/2) f Id^i is the mean intensity; S is 
the Lya emission function; a = H s v a jc is the increasing rate of the frequency redshift per 
unit distance in the z-direction (i.e. the direction in which the H I medium in a H I halo 
follows the pure Hubble expansion); and 

1>{e,w)= [£ 2 +(l-£> 2 ] 1/2 , (3) 

represents the effect of the anisotropic expansion of a H I halo where w is the z-axis compo- 
nent of direction cosine. The source function on the right-hand side of equation (2) can be 
written as 

N 

S = —^- 2 5(u)5(r), (4) 
(Air) r l 

where N a is the steady emission rate of Lya photons by the source (in photons s _1 ) and 
assumed to be constant. 

In the whole calculation in this paper, we assume that a steady point source galaxy 
at redshift 10 which lit up before the cosmological reionization is surrounded by a uniform, 
completely neutral IGM. In other words, we assume that the reionization epoch is redshift 
~ 10, which is intermediate between the late epoch of ~ 6 (Becker et al. 2001) and the early 
epoch of ~ 20 (Kogut et al. 2003). Apparently, the assumption of completely neutral IGM 
at redshift ~ 10 is oversimplified since a local cosmological H II region can spread over the 
source of Lya photons (e.g., Haiman 2002; Santos 2004). Although the so-called "proximity 
effect" is important for quasars but is less so for galaxies, the typical sizes of the H II regions 
around Lya source galaxies are expected to be not much less than ~ 1 proper Mpc, which 
is comparable to the size of the scattering region discussing in this paper (Barkana & Loeb 
2001). So there simply must be a phase when H II regions around early galaxies are becoming 
larger than the surrounding filaments; our results can hold prior to the phase. Fortunately, 
as shown in KK04, the Lya line profiles of the first galaxies at redshift > 10, which are 
considered to be in such phase, are analogous with those of sources at redshift 10 shown in 
this paper. However, some corrections for the model assumed in this paper due to the H II 
regions around Lja sources and the density fluctuations in H I halo may be necessary. 

Our calculation algorithm is identical with that of KK04. Avoiding redundancy, we only 
summarize the Monte-Carlo technique used in the calculations here (see §2 and §3 of KK04 
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for the details). At first, according to an initial frequency redshift of test photons from the 
Lya line center, the initial distances of them from the source are determined stochastically 
by means of the diffusion solution of Loeb & Rybicki (1999). The tentative optical depth r 
for a test photon is determined as — In i? sca , where i? sca is a random number between and 
1, and the photon proceeds the distance proportional to r. Once r > r max = J °° Xu+aipidl, 
we say that the photon escapes from the H I halo. During scattering events, the frequency 
redshift of a photon is increased according to the value of ip (e,w). These procedures are 
repeated until all test photons escape from the H I halo. 

The total number of the test photons N tot is 10 9 , which is ten times larger than that 
of KK04. Thanks to the increase of calculating photon number and the improvement of our 
numerical code, we can divide the whole solid angle into n regions with almost the same 
statistical precision as that of KK04. These regions are selected to have the same solid 
angle of An /n and to be axi-symmetric around the filament axis. Thus, it can be discussed 
meaningfully how much the Lya line profiles and the apparent Lya line luminosities of the 
first galaxies depend on the direction from which they are observed. We adopt z s = 10 and 
n = 10 for all calculations shown in this paper. 



3. RESULTS 

First, we define the Lya line profile per unit wavelength per unit solid angle (in sr _1 jimr 1 ), 
ip a (A, fi), as the following equation: 

Vai x,a)dXda = — . lhc/{1 + Zt)K] , (5) 

where N (A, fl) dXdVt is the number of escaping photons with the observed wavelength and 
solid angle in ranges of (A, A + dX) and (fl, Q + dfl), respectively; \ a = 0.1216 /iin is the res- 
onant wavelength of Lya in the rest-frame. The characteristic quantity on the denominator 
of the right-hand side of equation (5) represents the initial total energy of the test photons 
in the observer frame. Therefore, the integrated value of (p a (A, Q) over the whole wave- 
length and solid angle becomes unity only if the wavelength redshifts caused by scattering 
in expanding H I halos are negligible. 



3.1. Mean line profiles and scattered line luminosities 



Before we investigate the dependence of the observed Lya line profile and apparent Lya 
line luminosity of the first galaxies on viewing angles to them, we show the characteristic 
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effects of the anisotropic expansion of a H I halo to the Lya emission line according to KK04. 
For this purpose, here we define two tentative quantities: mean Lya line profile and scattered 
Lya line luminosity, denoted by y2™ ean (A) and I/™ ean , respectively. We emphasize here that 
they are not observational quantities except for those with e = 1 but make the effects of the 
anisotropic expansion of a H I halo clear. 

For a certain e, y2™ can (A) is obtained by counting the number of test photons with the 
observed wavelength in range of (A, A + dX) without distinction among escape directions, 
N(\)d\: 

{x) ~ N tot -[hc/(i + Zs) xy (6) 

Note that y?" lcan (A) is obtainable without knowing each line profile which is observed from 
a viewing angle. Figure 2 presents the mean Lya line profiles for the various values of the 
parameter e. The solid, long-dashed, short-dashed, dotted, long dash-dotted lines are e — 1.0 
(the Hubble expansion), 0.9, 0.7, 0.5, and 0.3, respectively. As e decreases, the wavelength 
at which y?™ can (A) has its peak value, A pea k, becomes longer as well as the peak value itself 
becomes lower; as a result, y2° lcan (A) becomes more diffuse in wavelength. These lines are 
well in agreement with the results of KK04. 

The scattered Lya line luminosity normalized by the total energy of the test photons 
in the observer frame, L™ ean , is obtained by the product of the integral of ip™ ean (A) over the 
whole wavelength and the whole solid angle 4n: 

L— = An <p~(\)d\. (7) 
Jxi 

This represents the total energy fraction of the scattered Lya line to the intrinsic one and 
therefore the value of L™ ean is always lower than (or perhaps equal to) unity. Recall that this 
is not an observational quantity for that with e ^ 1. A; and A u are selected to be 1.337 /im 
(= (1 + z s )\ a where z s = 10) and 2.00 jum, respectively; this wavelength range covers almost 
all photons escaped from the H I halo. The results of L" lcan are listed at the 7th column in 
table 1 for some e. One can see that L™ ean is always lower than unity and slightly decreases 
as e decreases. 



3.2. Characteristic line profiles and apparent line luminosities 

Here we present the numerical results of observed Lya line profiles and apparent Lya 
line luminosities of the first galaxies in order to investigate the dependence of them on 
viewing angles. First, we consider the observed line profiles of two characteristic cases; one 
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can be seen from z-direction and the other from r-direction denoted by ip z a (A) and <p r a (A), 
respectively. These are obtained by integrating (p a (A, Q) over the representative solid angle 
of An /n, and then, by multiplying the integrated value by a factor of nj 4n. We note here that, 
for e = 1.0, ip z a (A) and ip T a (A) completely correspond with each other and with ip™ ean (A). 
The results are depicted in figure 3 for e = 0.3 (thick and thin solid lines) and 0.5 (thick and 
thin dashed lines); the thick and thin lines are ip z a (A) and ip T a (A), respectively. y?™ ean (A) for 
e = 1.0 is also shown (dash-dotted line) for reference. One can see that, at the whole range of 
the observed wavelength, (p z a (A) with a certain e has always larger value than that of (p r a (A) 
with the same e. In addition, (p T a (A) rapidly fades down as e decreases, while cp z a (A) doesn't. 

The characteristic quantities of the observed Lya line profiles (i.e., A pca k and the FWHM) 
are summarized in table 1. At the 2nd and 3rd columns, A pea k of ip z a (A) and ip r a (A) with var- 
ious parameters of e are compiled, respectively. It is shown that A pca k of (p v a (A) is redshifted 
to longer wavelength than that of tp z a (A) with the same e. The 4th and 5th columns present 
the FWHMs of tp^ (A) and (p x a (A), respectively, and the 6th column shows the ratio of the 
FWHMs of tp z a (A) to (p r a (A). These indicate that (p r a (A) is always wider than ip z a (A) with the 
same e and that the difference among them grows further as e decreases. For convenience, 
we also give the measure of the observed flux per unit wavelength, f a (A), on the right-hand 
side of the vertical axis in figure 3 (in 10~ 20 ergs s" 1 cm~ 2 A^ 1 ), which is defined by 

msM) .-^-. (8) 

Here, di(z s ) is luminosity distance from the observer to the source redshift z s . For z s = 10 
which we adopt in this paper as the redshift of the Lya emission source, di is 3.27 x 10 29 
cm (~ 1.06 x 10 2 Gpc). We apply a typical Lya line luminosity of high-redshift LAEs of 
L a = 10 43 ergs s" 1 (e.g., Taniguchi et al. 2003; Hu et al. 2004; Malhotra & Rhoads 2004; 
Santos et al. 2004) as a reference luminosity of the Lya source at z s = 10; this is based on the 
fact that the luminosity function of LAEs in the redshift from 3 to 6 does not significantly 
change (Malhotra & Rhoads 2004; van Breukelen et al. 2005). The expected JWST limiting 
lOcr flux for a point source of 1 nJy to the source at redshift 10 with intrinsic Lya line 
luminosities of L a = 2.0 x 10 40 ,5.0 x 10 40 , 10 41 , 2.0 x 10 41 ,5.0 x 10 41 , 10 42 ergs s" 1 is also 
presented by contours in figure 3 (from top to bottom). 

The apparent Lya line luminosities can be obtained by integrating the observed line 
profile in a certain range of wavelength and by multiplying the integrated value by a factor 
of 4-7T like L™ ean in equation (7). This multiplication of the whole solid angle reflects that, 
when a Lya source at z s = 10 is observed, the Lya line luminosity is estimated under the 
assumption that the Lya emission from the source is isotropic with the observed flux density. 
Thus, this apparent Lya line luminosity coincides with the scattered Lya line luminosity 
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^mean Qn jy - m case that the Lya photons isotropically escape from the H I halo. In this 
paper, we adopt the integral interval from A; = 1.337 fim to \ u = 1.350 /im, assuming the 
source is observed with a tentative narrow-band filter centered at 1.344 fim with a bandwidth 
of 130 A. Then, we can obtain the apparent Lya line luminosities of L z a and L r a for ip^ (A) and 
tp T a (X), respectively. L z a and U a are listed at the 8th and 9th column in table 1, respectively, 
and their ratio is given at the last column. Compared with U a at the same e, L z a is always 
brighter. As e decreases, U a rapidly decreases while remains high value; therefore, the 
ratio of L z a to U a increases and results in the factor of > 43 for e < 0.3. 

3.3. Surface brightness distributions 

Our results also indicate that Lya emissions from the first galaxies are no longer point 
sources but diffuse on the sky; this agrees with the result of the previous work by Loeb & 
Rybicki (1999). This means that the surface brightnesses of the first galaxies at the Lya 
line become low. As well known, it is more difficult to detect more diffuse sources on the 
sky with same total luminosity (Yoshii 1993). In order to examine whether the predicted 
Lya source can be detected by the JWST, here we present the typical angular radius of the 
source at Lya line, in which the number fraction of incoming Lya photons becomes 0.90, in 
figure 4 for various e. The open circles with the thick solid line and open squares with the 
thick dashed line represent the typical angular radii of the source observed from the z- and 
r-directions, respectively. We also present the angular radius corresponding to the physical 
radius of O.lr* as the open triangles with the thin dotted line for reference. Numerical results 
are the open circles, open squares, and open triangles, while the various lines are the linear 
interpolators of them. Regardless of the observed direction, the typical angular size grows 
monotonically as e decreases; it ranges from 15" to > 50". The typical angular size of the 
source increases most rapidly in the case that the source is observed from r-direction. 

4. DISCUSSION 

According to the results described in the previous section, we discuss the effects of the 
anisotropic expansion of a H I halo to the observed Lya line profile and the apparent Lya line 
luminosity of a redshift 10 galaxy, and then, examine the detectability of the Lya emission 
by the JWST. 
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4.1. Detectability of Ly« emission: general comments 

As shown in figure 2 and summarized at the 7th column in table 1, the mean observed 
Lya line profile becomes more diffuse in wavelength and the scattered Lya line luminosity 
becomes dimmer as e decreases. This result implies that, on average, the detectability of the 
first galaxies by their Lya emissions gets worse as the expansion law of the H I halos becomes 
more anisotropic. Moreover, as e decreases, the discrepancy between the peak wavelength of 
the scattered Lya line and the resonant wavelength of Lya in the observer frame becomes 
more significant. That is, A pea k of the Lya line of the first galaxies does not correspond 
to (1 + z s ) \ a . This discrepancy causes a blunder of the Lya source redshift. In order to 
estimate the source redshift correctly, other lines of the same source such as Ha need to be 
detected. 

However, once the object with a very broad emission line which is very diffuse on the 
sky is detected, the effect of scattering in a H I halo can be easily recognized. This is 
because the FWHM of the Lya line is anomalously broad and corresponds to a Doppler 
width at the temperature of ~ 10 s K; this is not easily attainable by any other physical 
processes. Although the features of the Lya emissions from the first galaxies have not been 
fully understood yet, our results will present useful information because any intrinsic profiles 
can be concealed by this modified line profile. 

4.2. Dependence of detectability of Lya emission on observed direction 

Turning now to the dependence of detectability of Lya emissions from the first galaxies 
at redshift 10 on viewing angles, we find the following anisotropic effects of the H I halos. 
First, similar to the mean Lya line, the observed Lya line profile in each viewing angle be- 
comes more diffuse in wavelength as e decreases; this makes the apparent Lya line luminosity 
less luminous as e decreases. Second, however, the dimming and diffusing rates of the Lya 
line depend on observed direction; the apparent Lya line luminosity falls most sharply in 
r-direction, while it hardly decreases in the ^-direction as e decreases. 

The differences of the Lya line profile and of the apparent Lya line luminosity among 
viewing angles are understood by the dependence of the optical depth for Lya photon on 
propagative direction due to the anisotropic expansion of the H I halo. That is, because 
of the anisotropy in the expansion law of the H I halo, the optical depth is the largest 
for photons which propagate into the r-direction, while it is the smallest for those which 
travel to the z-direction. This means that escaping into r-direction is the most difficult; it 
needs to experience more scattering events and to get larger cumulative wavelength redshift 
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by scattering compared to that into z-direction. Therefore, we conclude that the direction 
along the axis of the underlying CDM filament is the most profitable to detect the Lja 
emission lines of the first galaxies. 

In any case, the observed Lja line profiles and the apparent Lja line luminosities of 
the first galaxies can be quite different from their intrinsic profiles and luminosities. This 
modified Lya line profile also needs to be fully understood in order to derive some physical 
quantities of the first galaxies, e.g., star formation rates (SFRs) and/or the escape fractions 
of the ionizing photons from their observed Lya line profiles (e.g., Miralda-Escude 1998; Cen 
k Haiman 2000; Madau & Rees 2000; Mesigner et al. 2004). 

The flux density represented by equation (8) gives a characteristic value as following: 

/ Q (A) = 9.3 x 10- 21 ( — ^ T ) ( ^ \ ergs s" 1 cm" 1 A" 1 . (9) 

This typical value of 9.3 x 10~ 21 ergs s _1 cm" 1 A -1 , which corresponds to 56 nJy and the 
AB magnitude of 27.0 mag at 1.34 /im, is about one order of magnitude brighter than the 
limiting lOu flux of JWST for a point source at the North Ecliptic Pole 1 . However, as shown 
in figure 4, our results indicate that Lya emissions from the first galaxies are diffuse on the 
sky; their typical angular sizes are a function of e and viewing angle, and range from 15" to 
~ 50". 

Evaluating S/N for the predicted Lya emission with an intrinsic Lya line luminosity of 
the source of L a = 10 43 ergs s" 1 by utilizing the JWST Mission Simulator (JMS) 2 , we find 
that the S/N for the source with e < 1.0 is always smaller than 10 over a > 100, 000 second 
exposure. Thus, it might be an observational challenge to detect the Lya emission against 
the very bright zodiacal light, even though it is at the North Ecliptic Pole and observed 
from the direction along the filament axis at the source. Furthermore, in the direction along 
the filament axis, there are expected to be clusters of galaxies. The presence of them at the 
foreground of the first galaxies might make the detection of the Lya emissions more difficult 
because it is very difficult to discriminate an object with low surface brightness against the 
one with higher surface brightness at the foreground of it. 

However, these clusters of galaxies can boost the detectability of the background Lya 
emission by gravitational lensing. According to the recent successful observations of the 
high-redshift LAEs by gravitational lensing (Kneib et al. 1996; Ellis et al. 2001; Kneib et 



1 For more details, see http://www.stsci.edu/jwst/ 

2 JMS can be used at http://www.stsci.edu/jwst/science/jms/index.html 
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al. 2004; Pello et al. 2004; Egami et al. 2005), high magnifications (~ 10 - 100 times) can 
be occurred. This high magnification is occurred only if the background source is in the 
so-called "critical regions". Location of the region is precisely known for well-understood 
clusters and depends on a redshift range of the source (Kneib et al. 1996; Ellis et al. 2001). 
By using the singular isothermal sphere (SIS) lens model (e.g., Turner et al. 1984) with the 
Einstein ring radius of 2' and Gaussian profiles as the surface brightness of the source, it is 
shown that the predicted Lja emission source can be magnified by a factor of > 10. In this 
case, the resultant S/N reaches > 10 over a > 100,000 second exposure. Thus, we conclude 
that detection of the predicted Lya source at the redshift z s ~ 10 is feasible over this decade 
with the JWST. 



4.3. Implications of detecting Lya emission from the first galaxies 

We insist that the Lya line profiles of the sources at redshift > 10 are analogous with 
those of a source at redshift 10 presented in this paper. They almost coincide with each other 
in dimensionless frequency redshift space as shown in figure 6 of KK04. Therefore, in order 
to evaluate the Lya line luminosities of the source at redshift > 10 and the typical angular 
radius of it at the Lya line, it is sufficient to know the characteristic frequency redshift and 
the characteristic proper distance at each redshift given in KK04. 

Furthermore, our results indicate that there is a correlation between the apparent Ly« 
line luminosity of the first galaxies and the underlying CDM structures at the initial con- 
traction stage. It is also shown that this correlation is further enhanced if the Lya source 
is observed by using a narrow-band filter targeting to a specific redshift as the recent obser- 
vations do. If one of the first galaxies at redshift ~ 10 is detected, this correlation can give 
us useful information that the underlying CDM filament at the extremely early Universe 
(cosmological age is younger than 0.5 Gyr) almost faces us, although degeneracy with the 
intrinsic Lya line luminosity of the galaxy remains. 

Possibly, however, Lya sources at z s > 10 are too faint to be detected as single Lya 
emission sources, but rather the assembly of them can be detected as a diffuse background 
source at near-infrared wavelength. This may have already been detected as the cosmic 
near-infrared background (CNIB), reported by independent groups (Wright & Reese 2000; 
Cambresy et al. 2001; Matsumoto et al. 2004) and cannot be explained by normal galaxy 
populations (Totani et al. 2001). Thus, considered with our results, the CNIB may be useful 
to map the large scale structures at the formation epoch of the galaxies in the early Universe. 



5. SUMMARY 



We investigated a Lya line transfer problem in the anisotropically expanding H I halos 
surrounding the first galaxies which are the sources of Lya photons. Using the improved 
numerical code, we solved the problem with 10 times larger number of test photons than 
that of KK04, and examined how much the detectability of the first galaxies at Lya lines 
depends on viewing angles to them. We found that the observed profiles and the apparent 
luminosities of the Lya emission lines, and the typical angular sizes of the first galaxies at 
Lya line strongly depend on the inclination angle to the axis of the underlying CDM filament. 
These physical quantities of the predicted Lya emission lines of the first galaxies are compiled 
for various e. The direction along the filament axis is found to be more profitable to detect 
the Lya emission than a direction perpendicular to the axis as the H I halo is more strongly 
decelerated against the cosmological Hubble expansion. In the perpendicular direction to 
the axis, because of the very broadened line profiles and the very diffuse surface brightness 
distributions, there seems to be little hope to detect the Lya emissions from the first galaxies 
with a typical Lya line luminosity of LAEs at redshift 3 — 6 of ~ 10 43 ergs s _1 . On the other 
hand, along the filament axis, there is expected to be clusters of galaxies as the sources of 
gravitational lensing. Based on our estimations by using a simple lens model and a simple 
surface brightness profile for the source, the high magnification with a factor of > 10 is 
expected to the predicted Lya emissions. We conclude that this high-magnification will 
allow us to detect the Lya emissions from the first galaxies in this decade with the S/N 
> 10 by JWST. Moreover, in order to know their physical quantities like the SFRs and the 
escape fractions of the ionizing radiation, it is essential to understand the relation between 
the predicted Lya emission line profile and the ratio of the apparent to the intrinsic Lya 
line luminosity. 
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Table 1. Summary of the physical quantities of the predicted Lya line of the first galaxies. 



Ap 0a ]j (/im) a FWHM (A) Apparent Line Luminosities 13 



1.0 1.339 1.339 52.7 52.7 1.0 0.98 0.70 0.70 1.0 

0.9 1.339 1.340 68.3 72.0 1.1 0.97 0.71 0.62 1.1 

0.7 1.340 1.342 84.8 115 1.4 0.96 0.71 0.42 1.7 

0.5 1.342 1.347 118 213 1.8 0.95 0.66 0.18 3.7 

0.3 1.346 1.363 190 506 2.7 0.92 0.46 0.01 43.6 



a The wavelength at which predicted Lya line profiles have their peak values. 

b The apparent Lya line luminosities obtained by integrating the observed Lya line profiles with wavelength in certain ranges and by 
multiplying the integrated value by a factor of 4tt. See the text for the integral intervals adopted in this paper. 

c The ratio of the FWHM of <p T a (A) to that of <p a (X). 
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Fig. 1. — Schematic illustration of one of the first galaxies and the underlying CDM filament 
surrounded by an anisotropically expanding H I halo with uniform density distribution, z- 
and r-axis (all directions perpendicular to z-axis) are also depicted. The thick arrows around 
the galaxy represent the expanding speed of H I medium assumed in this paper; along r- 
axis, the H I halo decelerated to some extent compared to that along z-axis (pure Hubble 
expansion) owing to the gravitational attraction force of the CDM filament. 
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Fig. 2. — Mean line profiles of the scattered Lya line of one of the first galaxies at z s = 10, 
^mean ag a f unc ti ori Q f t\y e observed wavelength A Q b s (in jttm) for various parameters of e. 
They show the line profiles of e = 1.0 (solid line; corresponds to the pure Hubble expansion), 
0.9 (long dashed line), 0.7 (short dashed line), 0.5 (dotted line), 0.3 (long dash-dotted line) 
as labeled. 
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Fig. 3. — Some of the predicted Lya line profiles of the first galaxies at z s = 10. Two 
characteristic cases are depicted; tp z a (A) and ip T a (A), which are the profiles of the Lya source 
observed from z- and r-directions, respectively. The two solid lines and two dashed lines are 
e = 0.3 and 0.5, respectively. The two thick lines with higher peak values in each kind of lines 
are <p z a (A), and the other two thin lines are ip z a (A). <^a ean (A) for e = 1.0 (dash- dotted line) is 
also presented for reference. The measure of observed flux densities of the Lya line, /a(A b s ), 
are also shown on the right-hand side of the vertical axis (in 10 -20 ergs s _1 cm" " 2 A" 1 ) with 
an intrinsic Lya line luminosity of L a = 10 43 ergs s _1 . Contours show the expected JWST 
limiting 10<r flux for a point source with intrinsic Lya line luminosity of L a = 2.0 x 10 40 , 5.0 x 
10 40 , 10 41 , 2.0 x 10 41 ,5.0 x 10 41 ,10 42 ergs s -1 (from top to bottom). The shadowed region 
represents the bandpass of a tentative narrow-band filter adopted in this paper (130 A 
centered at 1.344 /zm). 
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Fig. 4. — Typical angular radii of the first galaxies at Lja line on the sky (in arcsec) as 
functions of the parameter e. The open circles with thick solid line and the open squares 
with thick dashed line are the typical angular sizes of the sources observed from the z- and 
r-directions, respectively. The open triangles with thin dotted line are the typical angular 
sizes of a characteristic distance from the source of O.lr* for a reference. Note that the open 
circles, open squares, and open triangles are the numerical results, while the various lines 
are the linear interpolators of them. 



